Gluon mass through massless bound-state excitations. 
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Recent large- volume lattice simulations have established that, in the 
Landau gauge, the gluon propagator is infrared-finite. The most natural 
i-C 1 way to explain this observed finiteness is the generation of a nonpertur- 

bative, momentum-dependent gluon mass. Such a mass may be generated 
; gauge-invariantly by employing the Schwinger mechanism, whose main as- 

, sumption is the dynamical formation of massless bound-state excitations. 

• In this work we demonstrate that this key assumption is indeed realized 

by the QCD dynamics. Specifically, the Bethe-Salpeter equation describing 
the aforementioned massless excitations is derived and solved under certain 
approximations, and non-trivial solutions are obtained. 

PACS numbers: 12.38.Aw, 12.38.Lg, 14.70.Dj 

1. Introduction 

It is by now a well-established fact that large-volume lattice simula- 
tions in the Landau gauge yield a gluon propagator that reaches a finite 
non- vanishing value in the deep infrared [IH1]. Without a doubt, the most 
physical way of explaining this observed finiteness is to invoke the mecha- 
nism of dynamical gluon mass generation, first introduced in the seminal 
work of Cornwall [5] , and subsequently studied in a series of articles [6H8] . In 
this picture the fundamental Lagrangian of the Yang- Mills theory (or that of 
QCD) remains unaltered, and the generation of the gluon mass takes place 
dynamically, through the well-known Schwinger mechanism |9 w llj, without 
violating any of the underlying symmetries (for related contributions see, 
e.g., D2I). 

The main purpose of this presentation is to report on recent work [13] , 
where the Schwinger mechanism in quarkless QCD has been examined. 
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Fig. 1. Pictorial representation of the Schwinger mechanism in SDE. White (black) 
circles denote fully dressed propagators (vertices); a gray circle attached to the 
external legs indicates that they are background gluons. 



Specifically, the entire mechanism of gluon mass generation hinges on the 
appearance of massless poles inside the nonperturbative three-gluon vertex, 
which enters in the one-loop dressed diagram Fig. [1] of the Schwinger-Dyson 
equation (SDE) governing the gluon propagator. These poles correspond to 
the propagator of the scalar massless excitation, and interact with a pair of 
gluons through a very characteristic proper vertex, which, of course, must 
be non vanishing, or else the entire construction is invalidated. The way to 
establish the existence of this latter vertex is by finding non-trivial solutions 
to the homogeneous Bethe-Salpeter equation (BSE) that it satisfies. 

2. Gluon mass and the BS wave- function. 

In the formalism provided by the synthesis of the pinch technique (PT) 
with the background field method (BFM), known in the literature as the 
PT-BFM scheme [2], the Schwinger mechanism is integrated into the SDE 
of the gluon propagator, given in the Landau gauge byj 

&f u {q)=b ab ^M) ; V(«) = -iP^(q)A(q 2 ), (1) 

through the form of the three gluon vertex T afJiU , connecting one background 
with two quantum gluons. This is accomplished by adding a new non per- 
turbative piece to the three gluon vertex, which can be cast in the form of 
Fig. [H by setting 

K™ n (q,r,p) = I«\q)(^5 b ^B™ n {q,r,p). (2) 

In this expression, the non perturbative quantity 

Bftuiq, r,p) = Bxg^y + B 2 q^q u + B-^p^py + B 4 r^q v + B 5 r^p u , (3) 
1 The usual transverse projector is defined as -P M „(q) = g M „ — q^q^/q 2 . 
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is the effective vertex describing the interaction between the massless exci- 
tation and two giuons. B^ v is to be identified with the "bound-state wave 
function" (or "BS wave function") of the two-gluon bound-state, which, as 
we will see shortly, satisfies a homogeneous BSE. In addition, i/q 2 is the 
propagator of the scalar massless excitation. Finally, I a (q) is the (non- 
perturbative) transition amplitude, allowing the mixing between a (back- 
ground) gluon and the massless excitation. 

One can show |13j that, in the limit of zero momentum transfer q 2 = 0, 
the relevant quantity to consider is the form factor B\ appearing in Eq. (|3|). 
Furthermore, due to Bose symmetry with respect to the interchange \i -B- v 
and p <r-> r, we must have Bi(q,r,p) = —Bi(q,p,r), which implies in the 
aforementioned kinematical limit that -Bi(0, — p,p) = 0. So, when q — > 0, 
what survives is the derivative B[, which can be related to the derivative of 
the effective gluon mass through the exact all-order relation, 

[m 2 (p)]' = -I(0)B[(p). (4) 



3. The Bethe-Salpeter equation 

The existence of B[ is of paramount importance for the mass gener- 
ation mechanism; essentially, the question boils down to whether or not 
the dynamical formation of a massless bound-state excitation of the type 
postulated above is possible. As is well-known, in order to establish the 
existence of such a bound state one must (i) derive the appropriate BSE 
for the corresponding bound-state wave function, B^, (or, in this case, its 
derivative), and (ii) find non-trivial solutions for this integral equation. An 
approximate BSE for the bound-state wave function B^ v is given by the 
following expression [see Fig. [2] 

B arnn = f B fc^ P{k + g )A£(*)K™£ . (5) 
J k 

As explained in [13], Eq. ([5]) is obtained from the full BSE satisfied by the 
vertex T a ^ u . We will next approximate the four-gluon BS kernel /C by the 
lowest-order set of diagrams shown in Fig. [21 where the vertices are bare, 
while the internal gluon propagators are fully dressed. Going to Euclidean 
space, we define x = p 2 , y = k 2 , and z = (p + k) 2 ; then, following the 
procedure outlined in [13], the BSE becomes 



d0 sin 4 9 cos ( 



z + 10(x + y) + ^(x 2 + y 2 + lOxy) 



A(z) .(6) 
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Fig. 2. The BSE for the bound-state wave function f? M „ and the Feynman diagrams 
considered for the BS kernel. 



As a further simplification, we approximate the gluon propagator A(z) 
appearing in the BSE of ([6]) [but not the A 2 (y)] by its tree level value, that 
is, A(z) = 1/z. Then, the angular integration may be carried out exactly, 
yielding 



a s c A ( r x 2 y{o, 2 *y zy 2 



+ 



dyB[(y)A 2 (y)y 3 + 



25 x 3x 2 

4 y 4y2 



(7) 



4. Numerical analysis 

Next we discuss the numerical solutions for Eq. (|7|) for arbitrary values 
of x. Evidently, the main ingredient entering into its kernel is the nonper- 
turbative gluon propagator, A(g). In order to explore the sensitivity of the 
solutions on the details of A(g), we will employ three infrared-finite forms 
focusing on their differences in the intermediate and asymptotic regions of 
momenta, given by 

AfV) = Q 2 + m 2 , (8) 

A 2 -V) = V 2 

A 3 -V) = ^(, 2 ) + , 2 

The first expression ([8]) corresponds to a tree-level massive propagator, 
where m,Q is a hard mass treated as a free parameter. On the left panel 
of Fig. [3l the (blue) dotted curve represents A\ for rriQ = 376 MeV. The 
second model Q contains the renormalization-group logarithm next to the 
momentum q 2 , where p is an adjustable parameter varying in the range 
of p € [2,16]. Notice that the hard mass itiq appearing in the argu- 
ment of the perturbative logarithm acts as an infrared cutoff. The (black) 



1 + 



96vr 2 



In 



f + prriQ 



+ ml , 



(9) 



1 , lggAgg . 
1 + ^6^ ln 



+ P1 m2 g(<l 2 ) 
p? 



(10) 
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dashed line represents Eq. ([9]) when p = 16, mo = 376 MeV, and fi = 4.3 
GeV. Finally, Eq. (jlOp is simply a physically motivated fit for the gluon 
propagator determined by the large- volume lattice simulations of Ref. [2j, 
where m^(q 2 ) = m 4 /(q' 2 + ^m 2 ) is a a running mass. The values of the 
fitting parameters are m = 520 MeV, g\ = 5.68, p\ = 8.55 and, p2 = 1.91, 
and the (red) continuous line on the left panel of Fig. [3] represents this 
propagator. Notice that, in all three cases, we have fixed the value of 
A _1 (0) = ml « 0.14. 

Our main findings may be summarized as follows. 

(a) In Fig. [3j right panel, we show the solutions of Eq. ([7]) obtained 
using as input the three propagators shown on the left panel. For the simple 
massive propagator of Eq. (JSj) , a solution for B[ is found for a s = 1.48; in the 
case of A2 given by Eq. ([9]), a solution is obtained when a s = 0.667, while 
for the lattice propagator A3 of Eq. (|10p a non-trivial solution is found when 
a s = 0.492. 

(b) Note that, due to the fact that Eq. ([7]) is homogeneous and (effec- 
tively) linear, if B[ is a solution then the function cB[ is also a solution, for 
any real constant c. Therefore, the solutions shown on the right panel of 
Fig. [3] corresponds to a representative case of a family of possible solutions, 
where the constant c was chosen such that B[(0) = 1. 

(c) Another interesting feature of the solutions of Eq. ((JJ is the depen- 
dence of the observed peak on the support of the gluon propagator in the 
intermediate region of momenta. Specifically, an increase of the support 
of the gluon propagator in the approximate range (0.3-1) GeV results in a 
more pronounced peak in B[. 

(d) In addition, observe that due to the presence of the perturbative 
logarithm in the expression for A2 and A3, the corresponding solutions B[ 
fall off in the ultraviolet region much faster than those obtained using the 
simple Ai of Eq. ©. 

5. Conclusions 

In this presentation we have reported recent progress [13] on the study of 
the Schwinger mechanism in QCD, which endows gluons with a dynamical 
mass in a gauge invariant way. This mechanism relies on the existence of 
massless bound-state excitations, whose dynamical formation is controlled 
by a homogeneous BSE. As we have seen, under certain simplifying assump- 
tions, this equation admits non-trivial solutions, thus furnishing additional 
support in favor of the specific mass generation mechanism described in a 
series of earlier works [6HB]. In the future, this framework may serve as a 
starting point towards a consistent first-principle treatment of glueballs as 
bound states of gluons, at nonvanishing momenta. Specifically, one may en- 
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Fig. 3. The three models for the gluon propagator (left) and the corresponding 
solutions of the BS equation for B[(x) (right). 

visage a combined BS and SDE analysis, in the spirit of similar approaches 
in the case of color singlets qq bound states, employed in the literature for 
the dynamical description of mesons. 
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